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The magnitude of the CKM matrix elements V„b can be extracted from inclusive 
semileptonic B decays in a model independent way pioneered by Kolya Uraltsev 
and collaborators. I review here the present status and latest developments in 
this field. 


1. A semileptonic collaboration 


My continuing involvement with semileptonic B decays is mostly due to a fortuitous 
encounter with Kolya Uraltsev in 2002. A group of experimentalists of the DELPHI 
Collaboration, among whom Marco Battaglia and Achillc Stocchi, had embarked in 
an analysis of semileptonic moments and asked Kolya and me to help them out. He 
was the expert, I was a novice in the field and things have stayed that way for a long 
time thereafter. The joint paper that appeared later that year 1 contained one of 
the first fits to semileptonic data to extract \V c b\, the masses of the heavy quarks 
and some non-perturbative parameters, and it used Kolya’s proposal to avoid any 
I/m c expansion [2], The next step for us was to compute the moments with a cut 
on the lepton energy [3], as measured by Cleo and Babar (4} [5]. Impressed as I was 
by Kolya’s deep physical insight and enthusiasm, I was glad that he asked me to 
continue our collaboration. Kolya thought that a global fit should be performed by 
experimentalists, but as theoretical issues kept arising he tirelessly discussed with 
them every single detail; the BaBar fit [6], where the kinetic scheme analysis of 1 
was extended to the BaBar dataset, and the global fit of 7 owe very much to his 
determination. 

Kolya’s patience in explaining was unlimited and admirable: countless times I 
took advantage of it and learned from him. Our semileptonic collaboration later 
covered perturbative corrections [8], the extraction of V u b 9J, and a reassessment 
of the zero-recoil heavy quark sum rule 10 11 . It was during one of his visits 


to Turin that he suffered a first heart attack, but it did not take long before he 
was back to his usual dynamism. Working with Kolya was sometimes complicated, 
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but it was invariably rewarding. He was stubborn and we could passionately argue 
about a single point for hours. For him discussion, even heated discussion, was an 
essential part of doing physics. I will forever miss Kolya’s passionate love of physics 
and his total dedication to science. They were the marks of a noble soul, a kind 
and discreet friend. 

In the following I will review the present status of the inclusive B —> X c £v 
decays, the subject of most of my work with Kolya, who was a pioneer of the 
field. Semileptonic B decays allow for a precise determination of the magnitude of 
the CKM matrix elements V c b and V u b, which are in turn crucial ingredients in the 
analysis of CP violation in the quark sector and in the prediction of flavour- changing 
neutral current transitions. In the case of inclusive decays, the Operator Product 
Expansion (OPE) allows us to describe the relevant non-perturbative physics in 
terms of a finite number of non-perturbative parameters that can be extracted from 
experiment, while in the case of exclusive decays like B —> or B —> ir£H 

the form factors have to be computed by non-perturbative methods, e.g. on the 


lattice. Presently, the most precise determinations of V c b (the inclusive one 12 
and the one based onB-> D*£v at zero recoil and a lattice calculation of the form- 
factor 13 ) show a ~ 3er discrepancy that does not seem to admit a new physics 
explanation, as I will explain later on. A similar discrepancy between the inclusive 
and exclusive determinations occurs in the case of V u y, 14]. It is a pity that Kolya 
will not witness how things eventually settle. 


2. The framework 


Our understanding of inclusive semileptonic B decays is based on a simple idea: 
since inclusive decays sum over all possible hadronic final states, the quark in the 
final state hadronizes with unit probability and the transition amplitude is sensitive 
only to the long-distance dynamics of the initial B meson. Thanks to the large hi¬ 
erarchy between the typical energy release, of 0 (to{,), and the hadronic scale Aqcd, 
and to asymptotic freedom, any residual sensitivity to non-perturbative effects is 
suppressed by powers of Aqcd /mb- 

The OPE allows us to express the nonperturbative physics in terms of B meson 
matrix elements of local operators of dimension d > 5, while the Wilson coefficients 
can be expressed as a perturbative series in a s 15 19 . The OPE disentangles the 


physics associated with soft scales of order Aqcd (parameterized by the matrix 
elements of the local operators) from that associated with hard scales ~ TO&, which 
determine the Wilson coefficients. The total semileptonic width and the moments 
of the kinematic distributions are therefore double expansions in a s and Aqcd/t 71{,, 
with a leading term that is given by the free b quark decay. Quite importantly, 
the power corrections start at 0(Aq CD /m^) and are comparatively suppressed. At 
higher orders in the OPE, terms suppressed by powers of m c also appear, starting 


with 0 (Aq CD /toj x Aq c D / m c) 20 . For instance, the expansion for the total 
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semileptonic width is 
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free quark decay width, p = and A ew ~ 1.014 the leading electroweak 

correction. I have split the a 2 coefficient into a BLM piece proportional to (3q = 
11 — 2/3 rif and a remainder. The expansions for the moments have the same 
structure. 

The relevant parameters in the double series of Eq. |l]) are the heavy quark 
masses mb and to c , the strong coupling a s , and the B meson expectation values of 
local operators of dimension 5 and 6, denoted by /i 2 , Pq, p 3 Dl p 3 LS . As there are only 
two dimension five operators, two matrix elements appear at 0(l/m^): 

/ 4(/4 = ( 2 ) 


2 M b 


( 3 ) 


where n = —iD, D^ is the covariant derivative, b v {x) = e~ lmbV ' x b(x) is the b field 
deprived of its high-frequency modes, and G the gluon field tensor. The matrix 
element of the kinetic operator, p\, is naturally associated with the average kinetic 
energy of the b quark in the B meson, while that of the chromomagnetic operator, 
Pq, is related to the B*-B hyperfine mass splitting. They generally depend on 
a cutoff p = O(lGeV) chosen to separate soft and hard physics. The cutoff can 
be implemented in different ways. In the kinetic scheme 


21 22 


a Wilson cutoff 

on the gluon momentum is employed in the b quark rest frame: all soft gluon 
contributions are attributed to the expectation values of the higher dimensional 
operators, while hard gluons with momentum |fc| > p contribute to the perturbative 
corrections to the Wilson coefficients. Most current applications of the OPE involve 
0{l/m\) effects 23 as well, parameterized in terms of two additional parameters, 
generally denoted by p 3 D and p 3 LS 


22 


All of the OPE parameters describe universal 
properties of the B meson or of the quarks and are useful in several applications. 

The interesting quantities to be measured are the total rate and some global 
shape parameters, such as the mean and variance of the lepton energy spectrum or 
of the hadronic invariant mass distribution. As most experiments can detect the 
leptons only above a certain threshold in energy, the lepton energy moments are 
defined as 
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where Eg is the lepton energy in B —>• XJLv, r E e >E cat is the semileptonic width 
above the energy threshold E cut and dT/dEg is the differential semileptonic width 
as a function of Eg. The hadronic mass moments are 


(™x) 


1 

r Ee>E cu t 


[ ™ 2 x 

a ~>E C ut 


dT 

dni 2 x 


dm y . 


(5) 


Here, dT/dm 2 x is the differential width as a function of the squared mass of the 
hadronic system X. For both types of moments, n is the order of the moment. For 
n > 1, the moments can also be defined relative to (Eg) and (m 2 x ), respectively, in 
which case they are called central moments: 


^1 (-Ecut) — (Eg) E t >E aut , 
hi(E cut ) = (M x ) E e >E cut , 


hs(E C ut) = ((Eg - (Eg)f '*) El>Eout ; (6) 

h 2 , 3 (Ecut) = ((M 2 X (M 2 x ))^) Et >E mt . (7) 


Since the physical information that can be extracted from the first three linear 
moments is highly correlated, it is more convenient to study the central moments 
ti and hi, which correspond to the mean, variance, and asymmetry of the lepton 
energy and invariant mass distributions. 

The OPE cannot be expected to converge in regions of phase space where the 
momentum of the final hadronic state is O(Aqcd) and where perturbation theory 
has singularities. This is because what actually controls the expansion is not mb 
but the energy release, which is O(Aqcd) in those cases. The OPE is therefore valid 
only for sufficiently inclusive measurements and in general cannot describe differ¬ 
ential distributions. The lepton energy moments can be measured very precisely, 
while the hadronic mass central moments are directly sensitive to higher dimen¬ 
sional matrix elements such as /. i 2 and p 3 D . The leptonic and hadronic moments, 
which are independent of |F c b|, give us constraints on the quark masses and on 
the non-perturbative OPE matrix elements, which can then be used, together with 
additional information, in the total semileptonic width to extract | V c b |. 


3. Higher order effects 


The reliability of the inclusive method depends on our ability to control the higher 
order contributions in the double series and to constrain quark-hadron duality vio¬ 
lation, i.e. effects beyond the OPE, which we know to exist but expect to be rather 
suppressed in semileptonic decays. The calculation of higher order effects allows 
us to verify the convergence of the double series and to reduce and properly esti¬ 


mate the residual theoretical uncertainty. Duality violation, see 24 for a review. 


is related to the analytic continuation of the OPE to Minkowski space-time. It can 
be constrained a posteriori, considering how well the OPE predictions fit the ex¬ 
perimental data. This in turn depends on precise measurements and precise OPE 
predictions. As the experimental accuracy reached at the B factories is already 
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better than the theoretical accuracy for most of the measured moments and will 
further improve at Belle-II, efforts to improve the latter are strongly motivated. 

The main ingredients for an accurate analysis of the experimental data on the 
moments and the subsequent extraction of | V C 61 have been known for some time. 
Let us consider first the purely perturbative contributions. The 0(a s ) perturbative 
corrections to various kinematic distributions and to the rate have been computed 
long ago. In particular, the complete 0(a s ) and O(a 2 /3o) corrections to the charged 
leptonic spectrum have been first calculated in [25||27 and 28 . The so-called 
BLM corrections 29 , of O(a 2 /3o), are related to the running of the strong coupling 


inside the loops and are usually the dominant source of two-loop corrections in B 
decays. The first 0(a s ) calculations of the hadronic spectra appeared in 30 32 and 
were later completed in [§J 33; [34], while the O(a 2 f3o) contributions were studied 
in |8l [3l] 32] 34|. The triple differential distribution was first computed at 0(a s ) 


33 


its 0(a™/3g : ) corrections can be found in [8]. 


The complete two-loop perturbative corrections to the width and moments of 
the lepton energy and hadronic mass distributions have been computed in 


35- 


37 


by both numerical and analytic methods. The kinetic scheme implementation 

In general, using a s (mb) in the one- 


for actual observables can be found in 38 


loop result and adopting the on-shell scheme for the quark masses, the non-BLM 
corrections amount to about —20% of the two-loop BLM corrections and give small 
contributions to normalized moments. In the kinetic scheme with cutoff /r = lGeV, 
the perturbative expansion of the total width is 


T[B -> X c ev] oc 1 - 0.96 

+a82 (v) 


(m b ) 


-0.48 
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+ Q(a 3 s ) « 0.916 


( 8 ) 


39 


and 


Higher order BLM corrections of 0(o i ^I3q~ 1 ) to the width are also known ||8 
can be resummed in the kinetic scheme: the resummed BLM result is numerically 
very close to that of from NNLO calculations 39 . The residual perturbative error 
in the total width is about 1%. 

In the normalized leptonic moments the perturbative corrections cancel to a large 
extent, independently of the mass scheme, because hard gluon emission is compara¬ 
tively suppressed. This pattern of cancelations, crucial for a correct estimate of the 
theoretical uncertainties, is confirmed by the complete 0(a 2 ) calculation, although 
the numerical precision of the available results is not sufficient to improve the overall 
accuracy for the higher central leptonic moments 38 . The non-BLM corrections 


turn out to be more important for the hadronic moments. Even though it improves 
the overall theoretical uncertainty only moderately, the complete NNLO calculation 
leads to the meaningful inclusion of precise mass constraints in various perturbative 
schemes. 

The coefficients of the non-perturbative corrections of 0(Aq CD /to£) in the dou¬ 
ble series are Wilson coefficients of power-suppressed local operators and can be 
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Fig. 1. One-loop diagrams contributing to the current correlator. The background gluon can be 
attached wherever a cross is marked. 


computed perturbatively. The calculation of the 0(a s Aq CD /TO^) corrections has 
been recently completed. The 0(a s ) corrections to the coefficient of fi 2 have been 
computed numerically in 1401 and analytically in 41 . They can be also obtained 


from the parton level 0(a s ) result using reparameterization invariance (RI) rela- 
In fact, these RI relations have represented a useful check for the 


16 42 43 


tions 

calculation of the remaining 0(a s Aq CD /mj) corrections, those proportional to /i^., 
which was completed in [44| . The calculation consists in matching the one-loop di¬ 
agrams in Fig. [T] representing the correlator of two axial-vector currents computed 
in an expansion around the mass-shell of the b quark, onto local HQET operators. 

A recent independent calculation 45 of the semileptonic width at m c = 0 
seems to be in agreement with the m c —> 0 limit of 44 . Refs. 41 44 provide 
analytic results for the 0(a s Aq CD /mj) corrections to the three relevant structure 
functions and hence to the triple differential semileptonic B decay width. The most 
general moment have now been computed to this order and employed to improve 
the precision of the fits to \V c b\ fl2| . 

Numerically, using for the heavy quark on-shell masses the values mb = 4.6 GeV 
and m c = 1.15 GeV, the total semileptonic width reads 

,2 x ~ 


r B^XJv — To 


O- 178 ?) 
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2 ml 
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- 1.94 + 2.42 
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where To is the tree level width and we have omitted higher order terms of 0{a 2 s ) 
and 0{l/m\). The coefficient of /j, 2 is fixed by RI (or equivalently, by Lorentz 
invariance) at all orders. The parameter /Iq is renormalized at the scale mb- It is 
advisable to evaluate the QCD coupling constant at a scale lower than mj. If we 
adopt a s = 0.25 the 0(a s ) correction increases the Hq coefficient by about 7%. In 
the kinetic scheme with cutoff fx = lGeV and for the same values of the masses the 
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Fig. 2. Relative NLO correction to the pQ coefficients in the width (blue), first (red) and second 
central (yellow) leptonic moments as a function of the renormalization scale of p^,. 


width becomes 

r B^fX a iv = To 
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where the NLO corrections to the coefficients of are both close to 15% but 

have different signs. Overall, the 0(a s Aq GTl /ml) contributions decrease the total 
width by about 0.3%. However, NLO corrections also modify the coefficients of 
/rj, / Iq in the moments which are fitted to extract the non-perturbative parameters, 
and will ultimately shift the values of the OPE parameters to be employed in the 
width. Therefore, in order to quantify the eventual numerical impact of the new 
corrections on the semileptonic width and on | |, a new global fit has to be per¬ 

formed. The size of the 0(a s fj,Q/m%) corrections depends on the renormalization 
scale p, of the chromomagnetic operator. This is illustrated in Fig. [2] where the size 
of the NLO correction relative to the tree level results is shown for the width and 
the first two leptonic central moments at different values of /r. The NLO corrections 
are quite small for /i ss 2GeV and, as expected, increase with p,. For fj, <; the 
running of h'q appears to dominate the NLO corrections. In view of the importance 
of 0(l/m\) corrections, if a theoretical precision of 1% in the decay rate is to be 
reached, the 0(a s /m%) effects need to be calculated. 

As to the higher power corrections, the 0{l/m\) and 0(l/?ng) effects were 
computed in 46 . The main problem here is the proliferation of non-perturbative 
parameters: as many as nine new expectation values appear at 0(l/mf) and more 


at the next order. Because they cannot all be extracted from experiment, in 46 


they have been estimated in the ground state saturation approximation, thus re¬ 
ducing them to products of the known 0(l/m^’ 3 ) parameters, see also 


47 


In 


this approximation, the total 0(l/mg 5 ) correction to the width is about +1.3%. 
The 0(1/?71 q) effects are dominated by 0(l/?7i 3 m;?) intrinsic charm contributions, 


amounting to +0.7% 20 . The net effect on | F ch | also depends on the corrections to 
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the moments. Ref. 46 estimate that the overall effect on | V c b\ is a 0.4% increase. 


While this sets the scale of higher order power corrections, it is as yet unclear how 
much the result depends on the assumptions made for the expectation values. A 
new preliminary global fit [48] performed using different ansatz for the new non- 
perturbative parameters seems to confirm that these corrections lead to a small 
shift in | V c b \. 

Two implementations of the OPE calculation have been employed in global 
analyses; they are based either on the kinetic scheme [3j [2lJ [22, 39 or on the IS 1 
mass scheme for the b quark mass 49 50 . They both include power corrections 


up to and including 0(1/ml) and perturbative corrections of O(al/3o). Beside 
differing in the perturbative scheme adopted, the global fits may include a different 
choice of experimental data, employ specific assumptions, or estimate the theoretical 
uncertainties in different ways. Recently, the kinetic scheme implementation has 
been upgraded to include first the complete O(a^) 
contributions. 


38 


and later the (a s A 2 /m^) 


12 


4. | V c b | and the fit to semileptonic moments 

The OPE parameters can be constrained by various moments of the lepton energy 
and hadron mass distributions of B —> X c iv that have been measured with good 
accuracy at the B-factories, as well as at CLEO, DELPHI, CDF [6| p)l}|56|. The 
total semileptonic width can then be employed to extract \V c b\. The situation is less 
favorable in the case of \V u b \ 1 where the total rate is much more difficult to access 
experimentally because of the background from B —> X c £u, but the results of the 
semileptonic fits are crucial also in that case. This strategy has been rather success¬ 
ful and has allowed for a ~ 2% determination of V c b and for a ~ 5% determination 


of V u b from inclusive decays 14, 82 


The first few moments of the charged lepton energy spectrum in B —► X c £v 
decays are experimentally measured with high precision — better than 0.2% in the 
case of the first moment. At the B-factories a lower cut on the lepton energy, Eg > 
E cutl is applied to suppress the background. Experiments measure the moments at 
different values of E cut , which provides additional information as the cut dependence 
is also a function of the OPE parameters. The relevant quantities are therefore £ 1 , 2 , 3 , 
^■ 1 , 2 , 3 ) as well as the ratio R* between the rate with and without a cut 


R* (E cu t) = 


ft 


E c . 


dEt£r t 


r 

Jo 


dE 


( 9 ) 


This quantity is needed to relate the actual measurement of the rate with a cut to 
the total rate, from which one conventionally extracts \ V c f,\. All of these observables 
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Fig. 3. Two-dimensional projections of the fits performed with different assumptions for the theoretical 
correlations. The orange, magenta, blue, light blue 1-sigma regions correspond to the four scenarios 
considered in |[58j. The black contours show the same regions when the m c constraint of Ref. ([59] is 
employed. 


can be expressed as double expansions in a s and inverse powers of mb, schematically 




M™ + ( — Ym[ 2) +(M^’ 0) + 4 

V 7r / V tt / mr 


+ ( Mf’°) + A ^ + Mf' )P ^ + M^ P \ + ... (10) 

V n J m b m b m b 

where all the coefficients depend on m c , mb, E cut , and on various renormaliza¬ 
tion scales. The dots represent missing terms of O(af), 0(ot 2 s /m 2 ), 0{a s /m\), and 
0(l/m|), which are either unknown or not yet included in the latest analysis 12 . 
It is worth stressing that according to the adopted definition the OPE parameters 
/ij, ... are matrix elements of local operators evaluated in the physical B meson, 
i.e. without taking the infinite mass limit. 

The semileptonic moments are sensitive to a specific linear combination of m c 
and mb, « mb — 0.8m c 57 , see Fig. [3| which is close to the one needed for the 
extraction of | V c b \, but they cannot resolve the individual masses with good accu¬ 
racy. It is important to check the consistency of the constraints on m c and mb from 
semileptonic moments with precise determinations of these quark masses, as a step 
in the effort to improve our theoretical description of inclusive semileptonic decays. 
Moreover, the inclusion of these constraints in the semileptonic fits improves the 
accuracy of the | V u b\ and | V c b\ determinations. The heavy quark masses and the 
non-perturbative parameters obtained from the fits are also relevant for a precise 
calculation of other inclusive decay rates such as that of B —> X s y 
In the past, the first two moments of the photon energy in B —> 


58 


X s q have gen- 
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erally been employed to improve the accuracy of the fit. Indeed, the first moment 
corresponds to a determination of mb- However, in recent years rather precise de¬ 
terminations of the heavy quark masses (e + e _ sum rules, lattice QCD etc.) have 


become available, based on completely different methods, see e. g. 59-68 and 69 70 


for reviews. The charm mass determinations have a smaller absolute uncertainty 
and appear quite consistent with each other, providing a good external constraint 
for the semileptonic fits. Radiative moments remain interesting in their own respect, 
but they are not competitive with the charm mass determinations. Moreover, ex¬ 
periments place a lower cut on the photon energy, which introduces a sensitivity to 
the Fermi motion of the 6-quark inside the B meson and tends to disrupt the OPE. 
One can still resum the higher-order terms into a non-local distribution function 
and parameterize it assuming different functional forms 71 73 , but the parame¬ 
terization will depend on mb, etc., namely the same parameters one wants to 
extract. Another serious problem is that only the leading operator contributing to 
inclusive radiative decays can be described by an OPE. Therefore, radiative mo¬ 
ments are in principle subject to additional 0(AQCD/n r ib ) effects, which have not 


yet been estimated 74 . For all these reasons the most recent analyses 12 58 have 


relied solely on charm and possibly bottom mass determinations. 


The global fits of Refs. 12 58 are performed in the kinetic scheme with a cutoff 


p = 1 GeV and follow the implementation described in 3 38 . The two fits only 
differ in the inclusion of 0{a 2 s /m%) corrections and in the consequent reduction 
of theoretical uncertainties. In order to use the high precision m c determinations 
without introducing additional theoretical uncertainty due to the mass scheme con¬ 
version, it is convenient to employ the MS scheme for the charm mass, denoted by 
m c (p c ), and to choose a normalization scale p c well above m c , e.g. 3 GeV. 

The experimental data for the moments are fitted to the theoretical expressions 
in order to constrain the non-perturbative parameters and the heavy quark masses. 


43 measurements are included, see 58 for the list. The chromomagnetic expectation 
value Pq is also constrained by the hyperfine splitting 

1 


Mb * — Mb = - — + O 
3 mb 


a s p G 


m b m b , 

Unfortunately, little is known of the power corrections to the above relation and 


only a loose bound 75 can be set, see 11 for a recent discussion. For what concerns 
p 3 s , it is somewhat constrained by the heavy quark sum rules 
use the constraints 


75 


Refs. 


12 


58 


/4 = (0.35 ± 0.07) GeV 2 , p\ s = (-0.15 ± 0.10) GeV 3 


( 11 ) 


It should be stressed that p\ s plays a minor role in the fits because its coefficients 
are generally suppressed with respect to the other parameters. 

It is interesting to note that the fit without theoretical uncertainties is not good, 
with x 2 /dof ~ 2, corresponding to a very small p -value and driven by a strong 
tension (~ 3.5<r) between the constraints in Eq. (11) and the measured moments. 
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On the other hand, the fit without the constraints (11) is not too bad. Indeed 


theoretical uncertainties are not so much necessary for the OPE expressions to fit 


the moments — that would merely test Eq.(10) as a parameterization; they are 


instead needed to preserve the definition of the parameters as B expectation values 
of certain local operators, which in turn can be employed in the semileptonic widths 
and in other applications of the Heavy Quark Expansion. 

As noted above, the OPE description of semileptonic moments is subject to two 


sources of theoretical uncertainty: missing higher order terms in Eq. (10) and terms 


that violate quark-hadron duality. Only of the first kind of uncertainty is usually 
considered: the violation of local quark-hadron duality would manifest itself as an 
inconsistency of the fit, which as we will see is certainly not present at the current 
level of theoretical and experimental accuracy. 

In 12 we assume that missing perturbative corrections can affect the coeffi¬ 
cients of and /Mq at the level of ±7%, while missing perturbative and higher 
power corrections can effectively change the coefficients of p% and p\ s by ±30%. 
Moreover we assign an irreducible theoretical uncertainty of 8 MeV to the heavy 
quark masses, and vary a s (m&) by 0.018. The changes in M* due to these variations 
of the fundamental parameters are added in quadrature and provide a theoretical 
uncertainty 5M* h , to be subsequently added in quadrature with the experimental 
one, 8M^ xp . This method is consistent with the residual scale dependence observed 
at NNLO, and appears to be reliable: the NNLO corrections and the 0(l/m b ’ 5 ) 
(using ground state saturation as in 46]) have been found to be within the range 
of expectations based on the method in the original formulation of [3j. 

The correlation between theoretical errors assigned to different observables is 
much harder to estimate, but plays an important role in the semileptonic fits. Let 
us first consider moments computed at a fixed value of E cut : as long as one deals 
with central higher moments, there is no argument of principle supporting a correla¬ 
tion between two different moments, for instance t\ and hi- We also do not observe 
any clear pattern in the known corrections, and therefore regard the theoretical 
predictions for different central moments as completely uncorrelated. Let us now 
consider the calculation of a certain moment M, ; for two close values of E cut , say 
1 GeV and 1.1 GeV. Clearly, the OPE expansion for M l (l GeV) will be very similar 
to the one for Mj( 1.1 GeV), and we may expect this to be true at any order in a s 
and 1/m;,. The theoretical uncertainties we assign to M,( 1 GeV) and Mj(l.l GeV) 
will therefore be very close to each other and very highly correlated. The degree of 
correlation between the theory uncertainty of Mi(E\) and M,( E 2 ) can intuitively 
be expected to decrease as \E\ — E 2 I grows. Moreover, we know that higher power 
corrections are going to modify significantly the spectrum only close to the end¬ 
point. Indeed, one observes that the 0(l/m b ' 5 ) contributions are equal for all cuts 
below about 1.2 GeV (see Fig.2 of 46 ) and the same happens for the 0(a s p^/nril) 
corrections [40 . Therefore, the dominant sources of current theoretical uncertainty 
suggest very high correlations among the theoretical predictions of the moments for 












February 16, 2015 1:52 


World Scientific Review Volume - 9.75in x 6.5in 


memorial 


12 


Paolo Gambino 


Table 1. Results of the global fit with m c constraint in the default scenario 
of [12] . All parameters are in GeV at the appropriate power and all, except ra c , 
in the kinetic scheme at /x = 1 GeV. The last row gives the uncertainties. 



m c ( 3 GeV) 

A 

Pd 

As 

P 3 ls 

BR cfe (%) 

io 3 |f c6 | 

4.553 

0.987 

0.465 

0.170 

0.332 

-0.150 

10.65 

42.21 

0.020 

0.013 

0.068 

0.038 

0.062 

0.096 

0.16 

0.78 


cuts below roughly 1.2 GeV. 

Various assumptions on the theoretical correlations have been tried. A 100% 
correlation between a certain central moment computed at different values of E cut 
has been assumed e.g. in 0 )- This is too strong an assumption, which ends up 
distorting the fit because the dependence of M t on E cut , itself a function of the fit 
parameters, is then free of theoretical uncertainty. A fit performed in this way will 
underestimate the uncertainties. Another possibility has been proposed in Ref. 50 , 
with the theoretical correlation matrix equal to the experimental one. 

Four alternative approaches for the theoretical correlations are compared in ^58;. 
Fig-0 shows some of the results of the fits performed with the four options for the 
theoretical correlations. The fits include the two constraints of Eq. (11). In general, 
the results depend sensitively on the option adopted. In the case of the heavy quark 
masses, which are strongly correlated, we observe large errors differing significantly 
between the various options, although the central values are quite consistent. The 
results for the non-perturbative parameters depend even stronger on the option. 
The inclusion of precise mass constraints in the fit decreases the errors and neu¬ 
tralizes the ambiguity due to the ansatz for the theoretical correlations. It also 
allows us to check the consistency of the results with independent information. The 
effect of the inclusion of a precise charm mass constraint in the semileptonic fit is 
illustrated in Fig. [3] As expected, the uncertainty in the b mass becomes about 
20-25MeV in all scenarios, a marked improvement, also with respect to the pre¬ 
cision resulting from the use of radiative moments 


14 . The inclusion of the m r 


constraint indeed stabilized the fits with respect to the ansatz for the theory corre¬ 
lations. On the other hand, there is hardly any improvement in the final precision 
of the non-perturbative parameters and of \V c b\. 

Figs. [4] show two examples of leptonic and hadronic moments measurements 
compared with their theoretical prediction based on the results of the fit of 58 
with theory uncertainty. As anticipated, theory errors are generally larger than 
experimental ones. The situation is similar also for the fits in 

uses m c (3GeV) = 0.986(13) GeV 


12 


; the results are 


The default fit of 

shown in Table |T| where the bottom mass is expressed in the kinetic scheme, 
m fcm(i Q e y)_ Most available mb determinations, however, use the MS mass ) 

which is not well-suited to the description of semileptonic B decays as the calculation 
of the width and moments in terms of m.b(Wib ) involves large higher order corrections. 
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Fig. 4. Default fit predictions for I\ and /12 compared with measured values as functions of E cu t- 
The grey band is the theory prediction with total theory error. Filled symbols mean that the point was 
used in the fit. Open symbols are measurements that were not used in the fit. BaBar data are shown 
by circles, Belle by squares and other experiments (DELPHI, CDF, CLEO) by triangles. 


Alberti et al. (2014) 
Kuehn et al. (2009) 
Bodenstein et al. (2012) 
HPQCD (2014) 
Penin, Zerf (2014) 
Beneke et al. (2014) 
HPQCD,NRQCD (2014) 
Hoang et al. (2012) 
Lucha et al. (2013) 


Fig. 5. Comparison of different m^m^) determinations |l2] |60| |68 |76f|81] . The dashed line denotes 
the error before scheme conversion. 


t --r~ 


4.10 4.15 4.20 4.25 

m b(m b ) (GeV) 


Since the relation between the kinetic and the MS masses is known only to O(a^), 
the ensuing uncertainty is not negligible. It has been estimated to be about 30 


MeV 38 


mf n ( 1 GeV) - m b (m b ) = 0.37 ± 0.03 GeV, 


leading to a preferred value 


m b (m b ) = 4.183 ± 0.037 GeV, 


in good agreement with various recent mn b determinations 60 68 76^81 , as illus¬ 


trated in Fig(5] Of course, one can also include in the fit both m c and m b determi¬ 
nations, but because of the scheme translation error in m b the gain in accuracy is 


limited 12 58 . 
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As already noted, the semileptonic moments are highly sensitive to a linear 
combination of the heavy quark masses. When no external constraint is imposed 
on m Cj b, the semileptonic moments determine best a linear combination of the heavy 
quark masses which is close to their difference, 


m£ in (l GeV) - 0.85 m c (3 GeV) = 3.714 ± 0.018 GeV,. (12) 


The value of \V c t>\ is computed using 


\Vcb\ 


\V c b \ 2 BR c fo 

Tn r° p E 
1 B 1 B^X c tv 


(13) 


with tb = 1.579(5)ps 14j. Its theoretical error is computed combining in quadra¬ 
ture the parametric uncertainty that results from the fit, and an additional 1.3% 
theoretical error to take into account missing higher order corrections in the ex¬ 
pression for the semileptonic width. It turns out that using m c (2 GeV) rather than 
m c (2GeV) leads to a better converging expansion for the width, with smaller the¬ 
oretical error for \V c b\, about 1%. However, the value of \V c b\ extracted in this way, 
42.01(68) x 10~ 3 is compatible with that in Table 1. The same holds if the kinetic 
scheme is used also for m c , in which case we obtain 4% = 42.04(67) x 10 -3 . 

The fits are generally good, with x 2 /d.o.f. ~ 0.4 for the default fit. The low 
X 2 of the default fit is due to the large theoretical uncertainties we have assumed. 
It may be tempting to interpret it as evidence that the theoretical errors have 
been overestimated. However, higher order corrections may effectively shift the 
parameters of the 0(l/m 2 ) and 0(1 /to 3 ) contributions. If we want to maintain 
the formal definition of these parameters, and to be able to use them elsewhere, we 
therefore have to take into account the potential shift they may experience because 
of higher order effects. 

The fits with a constraint on m c are quite stable with respect to a change of 
inputs. In particular, small differences are found when experimental data at high 
E cut are excluded, and when only hadronic or leptonic moments are considered. One 
may also wonder whether the inclusion of moments measured at different values of 
E cu t really benefits the final accuracy. It turns out that the benefit is minor but 
non-negligible. 

In the kinetic scheme the inequalities p 2 (p) > Pq{p), p%(p)> hold a t 

arbitrary values of the cutoff p. The central values of the fit satisfy the inequalities. 
The dependence of the results on the scales of a s and on the kinetic cutoff has 
Changing the scale of a s from mb to 2 GeV the value of |I4&| 


been studied in 12 


increases by only 0.5%, well within its error, while mb increases by less than 0.4%. 
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5. Conclusions 


We have seen that the most recent value of Kb extracted from an analysis of 
inclusive semileptonic B decays is 12] 


\V cb \ = (42.21 ±0.78) x 10 


-3 


(14) 


A competitive determination of Kb derives from a comparison of the extrapolation 
of the B —» D*lv rate to the zero-recoil point with an unquenched lattice QCD 
calculation of the zero recoil form factor by the Fermilab-MILC collaboration 13], 


|Kb| = (39.04 ± 0A9 exp ± 0.53/ at ± 0.19 QE d) x 10 


-3 


(15) 


where the error is split into experimental, lattice and QED components. The values 
in Eqs. ( 14|15 ) disagree by 2.9cr. This is a long-standing tension, which has become 
stronger with recent improvements in the OPE and lattice calculations. A few 
comments are in order: 

• The zero-recoil form factor of B —> D*lv can also be estimated using heavy 


quark sum rules, see 10 11 for a recent reanalysis. Although this method 


is subject to larger uncertainties and it is difficult to improve its accuracy, 
it leads to a \V c b\ compatible with the inclusive determination, |Kb| = 
40.93(1.11) x 10 -3 . There exist also less precise determinations of |Kb| 
based on the decay B —> Dlv, but they do not help resolving the issue at 
the moment, see 82 for a review. 

The extrapolation of the B —> D*lv experimental data to the zero-recoil 
point is performed by the experimental collaborations using the Caprini- 


Lellouch-Neubert parameterization 83 , based on HQET at next-to-leading 


order and expected to reproduce the form factor within 2% (not included 
in the present error budget). While this rigid parameterization with only 
two free parameters fits well the experimental data at w ^ 1, at the present 
level of precision its use to extrapolate the rate to zero recoil is questionable. 
Lattice calculations of the form factors at non-zero recoil are currently under 
way; they would allow us to avoid the extrapolation. 

It is also possible to determine Kb indirectly, using the CKM unitarity 
relations together with CP violation and flavor data, without the above 
direct information: SM analyses of this kind by the UTfit and CKMFitter 

~ and (41.4t 3 'j) x 10" 3 ~ 


-3 


84 


collaborations give (42.05 ± 0.65) x 10 
both closer to the inclusive value of Eq. (|T41. 


'85 


In principle, the discrepancy between the values of | Kb | extracted from inclusive 
decays and from B —»• D*lv could be ascribed to physics beyond the SM, as the 
B —> D* transition is sensitive only to the axial-vector component of the V — A 
charged weak current. However, the new physics effect should be sizable (8%), and 
it would require new interactions which seem ruled out by electroweak constraints 
on the effective Zbb vertex 86 . The most likely explanation of the discrepancy 
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between Eqs. ( 14|15 1 is therefore a problem in the theoretical and/or experimental 
analyses of semileptonic decays. 

I do not have space here to discuss the closely related determination of \V u b\ 
from inclusive semileptonic B decays without charm, which shows a similar puz¬ 
zling tension between inclusive and exclusive determinations. However, I cannot 
avoid stressing the central role played by Kolya from the beginning also in this 
field 9, 72jj87,[88 . The interested reader is referred to 


82 for a recent review. 
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